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Q ■ abstract 

■ ^He atoms dissolved in superfluid ^He may form dimers ('^He)2 in two-dimensional (2D) geometries. 

We study dimer formation in films of dilute ^He-^He mixture. After designing a schematic ^He-^He 
interaction potential we calculate the dimer binding energy for various substrates. It is shown that 
^He impurity states localized near the substrate give rise to the largest magnitudes of the binding 
energies. 
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1 Introduction 



^He atoms dissolved in bulk ^He have always been considered an ideal system for testing the Landau 
approach describing the macroscopic properties of quantum fluids in terms of elementary excitations 
(quasiparticles) At low enough temperatures one can neglect the contribution of the excitations 
pertaining to superfluid ^He (phonons and rotons). The system is then described in terms of an in- 
teracting gas of ^He quasiparticles imbedded in the uniform background of superfluid ^He. A bare 
quasiparticle has a particle-like energy spectrum ||2[ with an effective mass larger than the bare mass. 
The quasiparticle interaction at long distances is the bare ^He-^He interatomic interaction renor- 
malized by a phonon-induced term j^. At small distances a strong repulsive core is predominant. 
However, the overall effect turns out to be attractive, as was demonstrated on the basis of theoretical 
calculations |Q, |5| and experimental data ^ on the s-wave scattering length a (a < 0). 

Despite the fact that the effective attraction between two '^He quasiparticles is too weak to lead 
to a bound state in the bulk, such a bound state - a (^He)2 dimer - should certainly exist in helium 
systems with reduced dimensionality § |9[. It is well known that in one and two dimensions any 
attractive potential satisfying the perturbation theory criterion gives rise to a bound state (see e.g. 
Ref. [^]). The result does not apply, however, in the case of a nonperturbative potential like the ^He- 
^He interaction, with a strong repulsion at short distances. Nevertheless, one can rigorously prove 
that spinless (^He)2 dimer must exist in quasi-2D and quasi-lD geometries (like films and narrow 
capillaries of ^He) provided a is negative and the characteristic scale of confinement of the ^He atoms 
is much larger than the interaction range - which is of the order of \a\ - (see Q and appendix A). When 
the system is cooled down to temperatures lower than the corresponding binding energy, the single 
impurity quasiparticles form (^He)2 Bose-dimers, and the Fermi component of ^He is replaced by a 
new Bose quantum fluid of (^He)2 of reduced dimensionality. This phenomenon should strongly affect 
the macroscopic properties of dilute mixtures resulting in an extra superfluid transition, new features 
of the phase diagram, anomalous sound absorption, etc. This is the reason why calculating the 
binding energy is the main issue in the theory of dimerized ^He-^He solutions. To be specific we will 
concentrate on dilute ^He-^He mixtures in various 2D geometries {i.e. in films on various substrates). 

In ^He films, ^He impurities are localized in the direction normal to the substrate over a character- 
istic width w which can be much smaller than the film thickness. When extrapolated to small w (of 
the order of a few atomic layers) the theory of Ref. |8|, |9| provides an estimate for the binding energy 
which seems to be quite attainable for experimental studies. However, it is difficult to assess the 
accuracy of such an estimate because the theory in question certainly does not hold in such restricted 
geometry, i.e. in a real 2D situation (when w ~ o). In a pure 2D case one cannot a priori be sure if a 
bound state of two '^He atoms exists at all. Therefore, calculating the binding energy for various film 
thickness and different substrates may provide direct clue on which experimental conditions are best 
suited for the detection of (^He)2 dimers. Such calculations are the primary goal of this paper. Here 
we focus on two cases where ^He impurities adopt a 2D configuration : 

(i) The Andreev states of ^He at the free surface of ^He (see |11, 12, ^] and references therein). 

(ii) The ^He states localized at the interface between a ^He film and the substrate, as predicted in 

MM- 

Case (ii) seems particularly favorable for the creation of dimers because ^He quasiparticles near 
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the substrate have a large effective mass (due to the higher ^He density) and a narrow wave-function 
(localized within the first ^He liquid layer). The localization length (the width of the wave function in 
the direction across the film) is of the order of a few angstroms. In Sec. 2 some possible realizations 
of a quasi-2D ensemble of ^He impurities in superfluid ^He films on various substrates are reviewed. 
In Sec. 3 we derive a schematic ^He-^He interaction and then evaluate in Sec. 4 the corresponding 
(^He)2 binding energy. We give our concluding remarks in Sec. 5. 



2 ^He impurities in ^He films. 

We assume in the following that inhomogeneities of the substrate in the x — y plane parallel to the 
surface play a negligible role and we impose translational invariance on the helium densities (^He and 
^He). This does not imply that possible corrugation of the surface is neglected. Indeed, in that case, 
'^He atoms will certainly fill the dips of the surface until it is more or less flat ; the effective substrate 
to be considered is then a mixture of the original substrate and these helium atoms, the effect of which 
is to weaken the potential. 

The first question to address is to determine the state of liquid helium on a given substrate. The 
different situations to expect, as a function of increasing attractive strength of the substrate potential, 
are the following : non wetting, wetting with prewetting [16|, and solidification of one or two layers. 



We limit our study to cases where no solid layer forms near the substrate : there is a clear qualitative 
change of '^He impurity states when ^He solidifies ; the ^He atoms occupy substitutional states in 
solid ^He and move through the lattice via tunneling processes (see e.g. Ref. fl^). Our continuous 
description of the solid misses this phenomenon, since the coupling of ^He atoms with the ^He matrix is 
represented here through a density dependent ^He effective mass fitted on properties of liquid mixtures 
(see details in Ref. Q). 

Hence, in order to specify the validity of our approach we need a criterion assigning a solid or liquid 
character to each of the first ^He layers. This has been done in Ref. [18| where different substrates 



were characterized by their interaction with helium through two coefficients C3 and D 

Here D is the depth of the attractive part of Vgubi^) and C3 characterizes the van der Waals tail. 
According to Ref. one can draw in the — D plane the lines corresponding to solidification of 
the first and second layer (see Fig. 1). For completeness, we have also reproduced, from Ref. |p!6|, t he 
line separating the non-wetting region from the wetting one. Notice that the criterion used in is 
only approximate and also that there are large uncertainties in the parameters of the potentials (up 
to ~ 30 % on D ||l^)- As an extreme case, for hydrogen the value C3 = 1000 KA^ extracted from 



different experiments (see [20] and references therein) is significantly larger than the theoretical value 
of 360 KA^ from Ref. ||T^. So, the predictions are uncertain for cases close to a line : it may well be 
that the first layer does solidify on Mg, and that the second layer solidifies on MgO, Cu and Ag. 

We now come to the ^He impurity states. Calculations were performed, as in Ref |14], in the 
limit of one '^He atom. When considering finite ^He coverage, each state generates a 2D Fermi disc, 
characterized by a 2D Fermi momentum. Also, through self-consistency, the ^He profile, the ^He 
average field - and consequently the single particle states - depend on the ^He coverage (see Ref. 
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[^] for a study of finite '^He coverage on the ^He bulk surface). We shall neglect these effects in the 
following discussion, which is valid for small ^He coverage only. 

Variation of the density functional with respect to the impurity wave function (p leads to the 
following equation 



d ( 



dz \ 2m* {z) dz 



+ Ue.tiz)(l)iz)=e<j)iz) , (2) 



Uext{z) is the ^He mean field, see Ref. [14|. It comprises a term due to ^He-'^He interaction plus 
the substrate potential Vgub. The density-dependent effective mass is parametrized as : 



(3) 



2m* (r) 2m3 V Pac 

where = 0.062 and P4(r) is the local ^He density averaged over a sphere of radius 2.38 A (see 
Ref. [0]). This parametrization is fitted to the pressure dependence of the ^He effective mass in bulk 
hquid ^He [0. 

Then the effective mass M* of the ^He impurity in a given state is defined by : 

Jl_ = / dz—^^<l)\z) (4) 
2M* J-oo 2m* {zY ^ ' ^ ' 

One can check that an energy close to that obtained through Eq. (|2|) can be recovered from the 
effective hamiltonian containing the constant effective mass M* 

^ = 2AF + ^'^'^'^ ' 

where Uext is the effective ^He potential appearing already in Eq. (^). 

Fig. 2 shows the results for the ^He density profile and the '^He wave function near a Cs and 
a Li surface. Numerical results for other substrates are displayed in Table 1. The existence of 2D 
^He states near a weakly binding surface, although somewhat counter-intuitive, appears as a general 
feature of these interfaces. 

The physics behind this feature is best understood by considering a Lekner approach to the prob- 
lem. In this variational method, the ansatz for the wave function ^' of — 1 ^He atoms and one ^He 
atom is taken as 

*(l,2,...,iV) =0(l)M/o(l,2,...,iV) (6) 

where cj) denotes the wave function of the impurity and ^'o the ground state wave function of N 
■^He atoms. Variation of the average energy of the system with respect to (f) results in a one-body 
Schrodinger equation (without effective mass effects) 

- ^ + UU^) 4>{z) = e ct>{z) (7) 
where U^^i{z) is an external field due to both the substrate and the ^He environment, 

^-(^) = silly + - + " + 
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Here (^^{z) = y/p^iz), 114 is the ^He chemical potential, t^^z) the kinetic energy density in the ^He 
ground state and Vsubi^) is the substrate potential. 

In the case of a free surface originally considered by Lekner, Vgub = 0; then Eq. (|8|) provides a 
mechanism for the formation of the pocket of potential at the surface : the kinetic energy density 
T4 goes to zero faster than the term in cp'^. For a numerical study, one assumes that T4(z) can be 
parametrized in term of the local density P4:{z). The simple following form was proposed in Ref. 

©a 

(fife))" p) 

where the value tq = 13.3 K is chosen so that the binding energy of a ^He atom in the bulk ^He is 
—2.8 K; Po is the saturation density of liquid ^He and n = 1.76 has been related to the excess volume 



parameter and the compressibility in dilute ^He-^He mixtures |22|. 

The effective potential U^^j- resulting from the Lekner theory for a film of 0.6 A^^ on a Mg substrate 
in shown in Fig. 3. The average field Uext obtained using the density functional theory is shown for 
comparison. One sees that the Lekner field is slightly more attractive at the surface; the corresponding 
energy of the Andreev state is found to be —5.24 K (using the bare '^He mass) whereas the density 
functional result is —5.27 K, with M* = 1.34 7713 (the experimental result being e = (—5.02 it 0.03) 
K and M* /m^ = 1.45 ± 0.1, see Ref. [|^) . Close to the substrate the average fields given by 
the two theories show similar qualitative behaviour, with attractive wells in correspondence with the 
oscillations of the '^He density. The Lekner field is more attractive than the density functional one, 
and leads to e = —14.1 K (using the bare mass), while the density functional theory gives e = —5 K, 
with M* = 2.9 m3. We believe that the density functional results are more reliable, since the model 
is fitted to a larger number of bulk mixture properties. The case of a semi-infinite liquid in contact 
with a Cs substrate is shown in Fig. 4. The Lekner value is e = —6.3 K, compared to the density 
functional value e = —4.74 K, with M* = 1.51 7713. 

Note that the Lekner result could be improved by a modification of the parametrisation (^). Here 
we only want to point out that substrate states are not specific to the Density Functional method, 
but emerge from the general features of the substrate-helium interface. 

To summarize, we stress here is that the mechanism by which an Andreev state is generated at 
the free surface of liquid helium is also operating at the liquid substrate interface. On a free surface, 
a ^He impurity is bound by 5 K (Andreev state), i.e. by 2.2 K more than in the bulk. Now a weak 
binding substrate perturbs the energetics of the free surface in two competitive ways, namely : i) the 
wall produces a readjustment of the ^He density profile, which reduces the width of the ^He wave 
function, and this tends to increase the energy of the state (to make it less bound); ii) however the 
attractive substrate potential acts on the ^He atom also and tends to decrease the energy of the ^He 
state. Clearly, if the perturbation is small enough, the gap of 2.2 K between the unperturbed surface 
energy and the bulk energy will not vanish, therefore a bound state has to remain at the liquid-wall 
interface. 

That the argument remains valid also in the case of Mg, which produces quite a marked layering 
of the fiuid, may seem surprising. Whether we reach in that case the limit of the model remains to be 
seen, by comparing with more microscopic calculations or with experiment. 



As discussed in [|14[ the existence of a substrate state allows one understand some unexplained 



temperature dependence of third sound velocities reported in Ref. |25|. Several experimental tests 
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of its existence were proposed in Ref. [14, ITH]. Recently, experimental evidence was reported that 



■^He impurities have a bound state at the ^He liquid-solid interface |26|, with a binding energy in fair 
agreement with the prediction of Ref. In the following we focus on one of the most exciting 

consequences : the possible formation of (^He)2 dimers with a sizeable binding energy. 

3 A schematic ^He— ^He interaction 

The next step in evaluating the dimer binding energy is to choose a sensible effective interaction 
between the ^He quasiparticles. The requirements are three-fold : (1) the long range attraction is 
the bare term reduced by a factor o? where a is the excess volume parameter in ^He (see |^ and 
below); (2) as mentioned in the introduction, the short distance term is repulsive and equal to the 
bare potential; (3) the effective potential should reproduce the s-wave scattering length a of ^He in 
^He (a ~ -0.97 A, see Ref. §). 

Let us start with the first requirement ; it will lead us to postulate a generic form for the effective 
potential. Our derivation of the energy of two ^He atoms imbedded in ^He is patterned on the approach 
of Ref. [^]. We consider two ^He atoms located at points ri and r2. Atom i (i = 1,2) occupies a 
volume Q3 centered at (we denote it as f^3[rj]) of which the ^He atoms are expelled (see Fig. 5). 
So, if the total volume of the sample is 0, the "^He's will occupy a volume Q' = Q — 0,3[ri] — 03[r2]. If 
Pi is the mean "^He density, then one can consider that each ^He atom occupies a volume = I//O4 
and the excess volume parameter is by definition 

The total potential energy of the system is 

^ = ^33 + Eu + ^44 , (11a) 

with obvious notations (for instance E34 is due to the '^He-^He interaction). E can be separated in a 
constant {i.e. position independent) term plus a term which is an effective ^He-^He potential : 

S = C^* + V;//(lri -ral) . (lib) 

The separation ( |llbD is made unambiguous by imposing that 14// goes to zero at large distance 
(when the two ^He impurities are far apart). 

Let us now study the contributions of £'33, -E34 and £'44 to the effective potential V^ff- 
If two helium atoms interact via the bare Lennard- Jones potential Vlj, -E33 is simply 



i^33 = V^Lj(|ri-r2|) . (12) 



where the explicit form of Vlj is 



Vuir) = As 



r J \r J 



(13) 



with e = 10.22 K and cr = 2.556 A. 
The cross term £34 is 
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Es, = [ VLj{\r, - r|)^ + / VLj{\r2 - r|)^ . (14) 
Jn' "4 Jn' "4 

In the r.h.s. of ( |l^ ) the two integrals are equal. Let us focus on the first one for instance. 
The integration domain can be separated into two contributions : 0,' = (0, — 0,s[ri]) — ^}3[r2]- The 
first subdomain describes the energy due to the introduction of atom 1 alone in the ^He matrix. It 
corresponds to the "dressing" of the bare ^He particle and plays no role in the quasiparticle interaction 
{i.e. it does not depend on |ri — and can be included in the constant term of Eq. ( |llb| )). The 
second subdomain describes the interaction of particle 1 with fictitious '^He atoms occupying J73[r2] 
and brings a contribution to Ve/f, since it depends of the respective position of the ^He impurities. 
Hence £"34 contributes to 14jj(|ri — r2|) with the following term (the factor 2 comes from the sum of 
the two equal integrals appearing in (0)) : 

E34^-2f VLj{\ri - r\)^ ^ -2^VLj{\ri - r2\) . (15) 

JQ3[r2] "4 "4 



In the last term of (|T5|) we have replaced the integral by an approximate form valid only if |ri — r2| 
is large compared with the characteristic radius of O3. If not, the dressing of the quasiparticles might 
affect one another (see below). 

The term £^44 is 

E..^Uf V^ilr-.n'^ff^ (16) 

^ J Jn'xQ' i'4 "4 

As above, the integration domain can be separated into several subdomains. The position depen- 
dent part of £44 gives a contribution to Vejf which represents the interaction of fictitious ^He atoms 
occupying the volumes ri3[ri] and f23[r2] 



^44-/ / VLA\r-r'\)^^^(^)\Lj{\r,-r2\). (17) 

J Jf73[ri]xn3[r2] "4 i'4 \i'4/ 

As in (^) the last term of Eq. (|l7|) is a long distance approximation. 

Gathering the contributions ( [l2| , p!5| , p!7| ) we see as stated in the beginning of this section that the 
long range quasi-particle interaction is equal to the bare interaction reduced by a factor . Essentially 
the same result was obtained in |^] using thermodynamical arguments in momentum space. In the 
terminology of Ref. Q the contribution of £33 +£34 would correspond to the direct part of the effective 
interaction and £44 would give the phonon-induced term. Of course, derivations based on excluded 
volume arguments such as in or as presented here (from |^^) are only valid for long wavelengths, 
i.e. for large distances between the ^He atoms. Working in real space has the advantage of providing 
a simple way to build an effective interaction sensible also at short distance : if the two atoms get 
closer one can mimic the interaction of the dressed particles by introducing a correlation term g{r) 
describing phenomenologically effects such as the disturbance of the '^He cloud around a ^He atom by 
the other particle. We choose for the correlation function the following form (r = |ri — r2|) : 



Kr) = exp{- ( ^ V} . (18) 
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Then the contribution of E34 to V^// is approximatively given by its long range approximation 
(the r.h.s. of (|l^)) multiphed by a factor g{r). Similarly the contribution of £^44 is multiplied by g'^{r), 
leading to a total effective interaction 



K//(r) = [!-(! + a)gir)Y VLj{r) . (19) 

The short range and long range behavior of the effective interaction (|^) follow the requirements 
(1) and (2) stated in the beginning of this section. Requirement (3) will be fulfilled by a correct choice 
of the free parameter Tc in (|T8|). The value Tc = 3.684 A gives the correct s-wave scattering length 
a = -0.97 A (see Ref. Q). The potential is shown on Fig. 6 where it is compared to the bare 
Lennard-Jones interaction. We have checked that working with a more realistic bare interaction such 
as the Aziz potential does not affect the qualitative picture presented below. Also the potential in 
films should be different from the bulk interaction : at the free surface for instance, ripplon — and not 
phonon — exchange should dominate the long range behaviour. Or equivalently, the excess volume 
parameter should be replaced by an "excess surface parameter" . We will not discuss this effect in the 
paper. 

In order to test the sensibility of the final results to the effective interaction we designed another 
potential V^ff having the required properties and roughly imitating at long range the potential derived 
by Owen using the hypernetted chain approximation ||29D. The potential was chosen to be : 



Veffir) = 4e 



I - I - C I - I cos I 27r ^ 



(20) 



This form obviously fulfills requirements (1) and (2) above. The parameters C and R are chosen 
so that Veff has a zero at r = 4 A as in Owen's results (this imposes R = 5.776 A) and that the 
s-wave scattering length has the correct value (this fixes C = 1.145). The corresponding 14// is 
plotted on figure 6. Note that the potential of Ref. (2^ is much deeper than Veff (it has a minimum 
at approximatively 19 K) and this would favour dimer's creation. On the other hand Owen uses no 
effective mass and this goes against binding. Hence we cannot directly compare our approach with 



the one of Ref. |29]; Owen's work is taken here only as an inspiration for designing a new potential in 
order to test the sensitivity of our results to the effective ^He-^He interaction. 



4 (^He)2 dimers 

The Hamiltonian describing two '^He quasiparticles in the film has the form: 

2 2 

^ = ^ + ^ + Veffilri - r2|) + Ue.t{zi) + Ue,t{z2) , (21) 

where Uext is the ^He mean field due to both the substrate and the ^He film (such as shown on Figs. 
2,3 and 4) and M* the ^He effective mass for the surface (Andreev) or substrate state (c/ Table 1). It 
should be pointed out that the ^He particles described by Eq. ( |2l| ) may, in general, be in two different 
localized states in the z-direction, i.e. may belong to two different 2D continua. In this case the 
effective masses of both quasiparticles in the Hamiltonian ( [2l| ) may be quite different. As explained 
earlier we restrict ourselves to considering the case of two ^He particles in the same state concerning 
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motion along the z-axis, since it provides the largest binding energy. In Eq. (^Tj) and in the following 
Veff is used as a generic notation valid for both potentials V^jf and K//- The description might in 
general be improved by phenomenologically introducing a '^He concentration-dependent effective mass 
M* within the local density functional approach to ^He-^He mixtures [13|. However the uncertainty 
in evaluating the dimer binding energy are such that this correcting term goes beyond the accuracy 
of the theory (see below). 

We make the following ansatz for the wave function of the two ^He atoms : 



^^(ri,r2) = ^e*^ll-^llx(r||)./.(zi)</.(z2) 



(22) 



Eq. (|22D describes two fermions with opposite spins. Ry is the center of mass coordinate in the 
{x,y) plane. For a dimer at rest (as we consider in the following) the corresponding momentum Ky 
is zero. The variable ry is the relative coordinate in the plane. We have separated the {x, y) and z 
direction. The functions (j){zi) and (j){z2) describe the motion along the z-axis. We assume (which 
is exact to the first order of perturbation theory) that they are of the type shown on Fig. 2, not 
being affected by the coupling between the two atoms. So the main point here is that the interaction 
between '^He atoms just slightly disturbs the motion in the z-direction but entirely changes the relative 
motion in the plane and leads to a bound state. In other words a perturbation theory can be applied 
for describing the normal motion only but the transverse dynamics of particles should be determined 
from the exact Schrodinger equation. In practice this means that one has to solve the 2D Schrodinger 
equation with a potential averaged over the unperturbed wave functions of ^He quasiparticles in the 
z coordinate. The function (j) is determined by the solution the equation : 



j1 

2M* 



ext \ 



where Pz is the z component of the momentum {z = z\ or Z'l) 



P? 



+ 



Pi 



2M* 2M 



2M* 



+ 



2M* 



+ 



II 



pir 



2M ^ 2/i 



(23) 



(24) 



with M and /i being respectively the total effective mass (2M*) and the reduced mass (M*/2). Then 
writing the Schrodinger equation for the entire wave function multiplying by (p* {zi)(j)* {Z2) and 
integrating over the variables z\ and Z2-, one easily finds : 



X(r||) = EdimXi^w) ' 



(25) 



where 



(V;//(r||)) = / (l)\zi)(l,\z2)Veff{vi-r2)dzidZ2 



(26) 



According to our perturbation scheme, all the z-dependence has been removed from Eq. (^51). The 
substrate and the ^He density play an indirect role, through the determination of {i.e. M*) and the 
z-wave function. 
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It appears that for all the cases we are interested in (the Andreev state or the substrate state) the 
z-wave function can be represented to a fairly good approximation by a simple gaussian (characterized 
by its half-width). In order to explore the sensitivity of the dimer binding energy to the two most 
relevant parameters, namely the effective M* and the width of the ^He wave function, rather than 
considering only the substrates listed in Table 1, we will give results for various values of M* and 
normalized gaussians for (j)'^{z) with various half width. 

In the limit of a large half-width, although common perturbation theory cannot be directly applied 
when looking for the solution of the Schrodinger equation, one can use the Fermi renormalization 
technique (see [^] and appendix A) to obtain the binding energy with a logarithmic accuracy: 

|E.,™l--g^e,p(-^^) , (27) 

where ro is a quantity of the order of the interaction potential range, is the effective mass of a 
^He quasiparticle in bulk ^He, and L is defined as: 

i = | \m\'dz. (28) 

Here L is of the order of the half-width, w, and it is a measure of the spatial extension of the 
z-wave function. The result above is valid only for large values of L/\a\ (see Appendix). One can 
see from Eq. (^) that in this limit the requirements of Sec. 3 are enough to determine the order of 
magnitude of the dimer binding energy. 

In the opposite limit of very localized states — such as those we are primarily interested in — our 
numerical results show a great sensibility of the binding energy to the details of the potential (the 
results for E^im as a function of the half-width for different effective masses are shown on Fig. 7 for V^// 
and on Fig. 8 for 14//). For instance, for M* /m^ ~ 3.1 and w ~ 0.5 A the binding energies estimated 
with Veff and Veff differ from each other by two orders of magnitude. Under these conditions an 
accurate quantitative prediction would be illusive. We just note that if Veff mimics accurately the 
exact quasiparticle interaction, it would be very difficult to observe the formation of (^He)2 bound 
states in experiment. On the other hand, if the effective potential 14// can be applied in the case in 
question, then for a magnesium substrate {M* /m^ = 2.9 and w = 0.9 A) the dimer binding energy 
would be Edim = —1-1 mK, which is a reachable temperature with modern experimental techniques. 
A problem might occur though because the first '^He layer on a Mg substrate could be solid (see the 
discussion in Sec. 2). The next candidate for dimer formation would then be hydrogen. Using the 
two potentials designed in Sec. 3 we could not get a reasonable dimer binding energy, mainly because 
M* in this case is not large enough. However, other pseudo-potentials might give a different result. 

Note also that the dimer binding energies are very sensitive to the exact value of the s-wave 
scattering length. This can be seen in the large L (or equivalently large w) limit from Eq. (^). We 
verified numerically that this is also true for small L : taking the outdated value a = —1.5 A we 
obtained (using 14//) a dimer binding energy on magnesium S^jm = — 6.3 mK. Hence experimental 
information on dimer binding could give valuable insight on the exact value of a. 

Although not quantitatively predictive, our study allows us nevertheless to draw a very clear 
qualitative picture. In highly compressed substrate layers the ^He effective mass is large. This reduces 
the kinetic energy and favours the creation of ('^He)2 bound states. The small spatial extension of the 
wave function in these layers is also favorable to the formation of dimers as can be seen from Eq. (|27| ) 
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and Fig. 7 and 8. From these figures one can also see that there is an optimal width of the order of 
0.5 A. It is interesting to note that the more attractive the substrate, the closer one approaches this 
value (see Table 1). The proper way to observe dimers would then be to choose a substrate for which 
the first ^He layer in close to be solid but still remains liquid. 



5 Concluding remarks 

After having classified various substrates according to their ability to solidify one or two layers of a 
multilayer helium film, we have concentrated on cases in which helium remains liquid. In addition to 
the already known Andreev states we have shown that ^He impurities were able to form a new 2D 
Fermi system near the substrate. After designing a sensible schematic interaction we have computed 
the binding energy of (^He)2 dimers in the limit of small ^He concentration. We found that for 
attractive substrates (such as magnesium or hydrogen) there was a reasonable hope to form dimers 
with sizeable binding energy. Note here that the binding energy of dimers in Andreev surface states 
is extremely small according to all our estimates, so the eventual observation of (^IIe)2 dimers would 
be a direct consequence of the existence of ^He substrate states. 

In this paper the quantitative calculation of the dimer binding energy was based on a semi- 
empirical effective interaction and the one-particle wave functions obtained within a density-functional 
approach. It would be very useful to carry out the appropriate calculations with the help of other 
theoretical approaches {e.g. Ref. [^]). This would allow a better understanding of the accuracy of 
the present computations. 

We hope that these results will motivate experimental study of the substrate ^He states as was 
proposed earlier in [0, The eventual formation of (^He)2 dimers would be a very interesting conse- 
quence of the existence of these states leading to an amazingly rich phase diagram. An experimentalist 
could face an extra Kosterlitz-Thouless phase transition, liquefaction of (^He)2 or polymerization of 
■^He, crossover from a Bose gas of dimers to a 2D Fermi fluid (with strong pair correlation) upon 
increasing the ^He concentration, etc. ^, 33, The dimer binding energy is a cornerstone 



characteristic of all these phenomena which are exciting objectives for experimental and theoretical 
studies. 

A natural continuation of the present work would be to study the binding of two dimensional 
clusters of ^He atoms. Although Fermi statistics favours droplets with even numbers of ^He atoms, 
there might be no upper limit for N , and our results could be a hint of the existence of a 2D liquid 
phase of ^He. On the basis of a regular arrangement in two dimensions (with 6 nearest neighbours 
per atom) we estimate the saturation energy of the liquid phase to be (in the most favourable case) 
of the order of 5 to 6 mK. This is to be related to the recent finding of Brami et al. ||3^ whose 
variational results lead to propose at low temperature a new "self-condensed" fluid phase for pure ^He 
films on graphite. As the presence of "^He favours the formation of ^He dimers it could also favour the 
formation of a liquid phase. 
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Appendix A 



In this appendix we derive the expression ( |27D for the dimer binding energy -E^im using the Fermi 
renormahzation technique (see |^^). According to this method we introduce a weak two-particle 
pseudo-potential Vf{r) which is assumed to meet the perturbation theory criterion. Let us emphasize 
that such an pseudo-potential has nothing to do with the interaction between quasiparticles which, 
indeed, cannot be treated in terms of a perturbation theory at all. 

The pseudo-potential Vf is supposed to result in the true scattering amplitude when calculated 
within the Born approximation : 



(Al) 



where is the effective mass of a ^ 



He quasiparticle in bulk ^He (M*^ = 2.887713). The idea of the 
method is to carry out all calculations in the framework of a perturbation theory for Vf and then to 
express the final formulae through the true s-wave scattering length a only by means of the renormal- 
ization ([All) . Thus the method works while all obtained expressions contain the pseudopotential Vj 
only in the integral form ([Al]). 

The solution of the Schrodinger equation, i.e. of Eq. ( p5| ) with Vj replacing Vg//, may be expressed 
in the form (see e.g. Ref. [fH, |3l[) 



Edi; 



TT exp 

2^r2 



{Vf)pdp 



-1 



(A2) 



where rg ~ |a| is of the order of the interaction range and (Vf) is defined by ( pq ) with Vf replacing 
Veff. Let us recall once again that the pseudopotential Vf is assumed to meet the criterion of a 
perturbation theory : for instance the integral entering Eq. 
integral as I one can write : 



([A2|) should converge. Denoting this 



{Vf)pdp = / Vf{^p^ + z^)F{z)pdpdz 



(A3) 



where 



F{z) = / {(Pit - z/2)|2 [0(t + z/2)fdt . 



(A4) 

Let the quantity w be a characteristic localization range in the z-direction for the wave function 
(j) {w describes the half-width introduced in section 4). If the localization length is sufficiently large, 
w ^ rQ, one can easily find that 



/ = ^ / Vf{r)d'r with i EE F(0) = J |0(t)|^di . 



(A5) 



Evidently L is of the same order as w (if 0^ is a normalized gaussian of half- width w, then 
L = i(;-y7r7ln4). Combining Eq. ( [A5| ) and the renormahzation relation ( [Al| ) yields : 



E, 



dim 



M*rl 



exp 



L 

M* |a| 



(A6) 
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In ( |A6D we have replaced the reduced mass by its value for the state (j){z) considered : fi = M*/2. 
We see that the binding energy is expressed in terms of the true scattering amplitude and does not 
contain a pseudopotential in any explicit form. Thus using the Fermi method is completely justified 
in the limiting case w ^ rQ. Note, however, that it is difficult to exactly determine the quantity ro 
(see 1 31]) which could, in general, be pseudo-potential dependent. On the other hand, its order of 



magnitude is known (ro ~ \a\) and, in fact, it does not enter the exponential factor which determines 
the dominant behavior of Edim- 

In the limiting case w <^ ro the integral / obviously reduces to 

I = jVf{p)pdp, (A7) 

and the binding energy explicitly depends on what kind of a pseudopotential is chosen. It means that 
the renormalization method is no longer valid and the 2D Schrodinger equation with the real interaction 
potential should be solved to find the binding energy. In this paper we have chosen to perform the 
calculation of E^nm by means of semi-empiric pseudopotentials which gave rise to reasonable estimates 
for quantities measured in experiment (see Sec. 2). 
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Table captions 



Table 1. Energy (e) , effective mass (M*) and half- width (w) of the ^He impurity state on various 
substrates and for the Andreev state. In the case of H2, the value C3 = 1000 KA^ extracted from 
different experiments |^0| is different from the theoretical value of 360 kA'^ [p^]. 



substrate 


C3 [ kA3 ] 


D[K] 


e [K] 


M*/m3 


w[A] 


Cs 


673 


4.4 


-4.74 


1.51 


2.0 


Rb 


754 


5.0 


-4.23 


1.78 


1.9 


K 


812 


6.3 


-4.15 


1.84 


1.8 


Na 


1070 


10.4 


-4.04 


2.02 


1.4 


Li 


1360 


17.1 


-4.17 


2.28 


1.2 


Mg 


1780 


32.0 


-5.00 


2.90 


0.9 


H2 


360 


28.0 


-5.51 


2.41 


0.9 




1000 


28.0 


-5.00 


2.61 


0.9 


^He free surface 






-5.20 


1.35 


3.6 
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Figure captions 



Figure 1. Phase diagram of ^He on various substrates at zero temperature (form [|T8|). The 



dotted hne is the wetting line |16] 



Figure 2. ^He hquid density (dotted Hne), ^He substrate state wave-function (with arbitrary 
units) and '^He average field on a Cs substrate (less attractive) and Li substrate (more attractive). 



Figure 3. ^He film (0.6 atoms per A^, i.e. 7.7 atomic layers) and ^He mean field on a magnesium 
substrate. The solid line is the Density Functional mean field Uext and the dashed line is the result of 
Lekner's approach U^^f.. 



Figure 4. Infinite ^He film and ^He mean field on a cesium substrate. The solid line is the 
Density Functional mean field Uext and the dashed line is the result of Lekner's approach U^xt- 



Figure 5. Schematic representation of the behaviour of two ^He atoms situated at points ri and 
r2 in bulk liquid ^He. Atom i {i = 1, 2) occupies a volume il3[rj] of which the ^He's are expelled (see 
the text). The shaded zone represents the volume occupied by liquid ^He. 



Figure 6. ^He interaction potentials. The solid line represents the bare Lennard-Jones potential 
Vlj, the dashed lines represent the effective potentials V^ff (long dashes) and V^ff (short dashes). 



Figure 7. Edim as a function of w computed with V^// for the effective masses M* /m^ = 3.1, 
3.15 and 3.2. 



Figure 8. Edim as a function of w computed with Vejf for the effective masses M* /m^ = 2.8, 
2.9 and 2.95. Note the change of energy scale with respect to Fig. 7. 
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